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FOREWORD 


This  worK  was  performed  for  and  funded  by  the  NAVSEA  Explosives  Development, 
Effects  and  Safety  block  program,  SF33-337-691.  The  results  and  conclusions 
presented  in  this  report  concerning  shock  waves  in  water  from  detonating  pentolite 
spheres  should  be  of  interest  to  those  seeking  underwater  sensitivity  information 
and/or  underwater  damage  information. 

The  authors  wish  to  acknowledge  the  assistance  of  J.  Marshall  and  A.  Brown 
in  conducting  the  lithium  niobate  gage  experiments;  D.  Gillmore  and  R.  Baker 
for  conducting  the  underwater  sensitivity  experiments  from  which  the  shock 
velocities  in  water  were  obtained.  H.  Jones  and  J.  W.  Forbes  obtained  the 
gage  information  on  the  6.84  cm  radius  pentolite  sphere.  Keith  Harrison  provided 
electronics  support  and  ran  many  of  the  computer  programs  used  for  data  reduction. 
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INTRODUCTION 


The  generation  of  spherical  shock  waves  in  water  is  of  considerable  interest 
to  researchers  and  theoreticians.  Spherical  shock-generating  systems  in  water 
have  been  used  to  study  the  behavior  of  shock  waves,  the  effect  of  the  shock 
waves  against  submerged  structures,  and  interactions  of  the  shock  waves  with 
other  spherical  shock  waves.  These  shock-generating  systems  also  have  been 
used  to  induce  detonation  in  explosive  samples  placed  at  various  distances 
from  a  spherical  donor. 2»3  such  a  system  was  used  by  Liddiard^  to  detect 
chemical  reaction  in  test  samples  at  pressures  as  low  as  4  kbar  in  the  sample 
or  3  kbar  in  water. 

In  applications  of  this  kind,  accurate  measurements  of  the  shock  pressure 
in  water  are  needed  for  reliable  assessment  of  results,  particularly  at  very 
close  range.  Hantel  and  Davis5  developed  and  cal Ibrated  a  spherical  shock¬ 
generating  system  using  carefully  manufactured  explosive  components,  the  main 
charge  being  composed  of  PBX  9205.  The  excellent  reproducibility  of  their 
system  would  make  it  ideal  for  use  at  very  close  distances  except  that  the 
cost  of  making  the  precision  charges  is  prohibitively  high  for  our  use  in 
explosive  testing.  {We  have  found  that  cast  pentolite  donor  charges,  made  with 
reasonable  care,  ^  te  very  good  results.) 

In  the  present  work,  the  shock  velocities,  pressures,  and  pulse  time 
histories  in  fresh  water  are  measured  as  a  function  of  distance  from  the  pentolite 

^Coleburn,  N.  L.  and  Roslund,  L.  A.,  “Interactions  of  Spherical  Shock  Waves  in 
Water,"  Proceedings »  Fifth  Symposium  on  Detonation,  18-19  Aug  1970,  pp.  581-588. 

Winning,  C.  H.,  "The  Underwater  Shock  Wave  Initiation  of  Cast  Pentolite,"  in 
Proceedings  of  the  Royal  Society,  Vol.  246,  No.  1245,  Jul  1958,  pp.  288-2%. 

talker,  F.  E. ,  and  Wasley,  R.  J.,  "Initiation  Patterns  Produced  in  Explosives  by 
Low-Pressure,  Long-Duration  Shock  Waves,"  Combustion  and  Flame,  Vol.  22,  1974, 
pp.  53-58. 

\iddiard,  T.  P.,  "The  Initiation  of  Burning  in  High  Explosives  by  Shock  Waves," 
Proceedings,  Fourth  Symposium  on  Detonation,  12-15  Oct  1965,  pp.  487-495. 

5Hantel,  L.  W.  and  Davis,  W.  C.,  "Spherical  Explosions  in  Water,"  Proceedings, 
Fifth  Symposium  on  Detonation,  18-21  Oct  1970,  pp.  599-604. 
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charge.  The  shock  velocity  as  a  function  of  distance  is  derived  from  framing 
and  smear-camera  records.  The  peak  pressures  and  pulse  time  histories  are 
measured  using  lithium  niobate  transducers.  The  pressure-vs-di stance  calibration, 
used  in  the  earlier  work  of  reference  4,  was  derived  from  shock  velocities 
obtained  from  a  relatively  few  framing-camera  records.  Since  that  work  was  done, 
many  more  data  have  been  obtained  and  refinements  have  been  made  in  analytical 
techniques. 

As  to  pressures  much  farther  away  from  the  charge,  tourmaline-gage  measure¬ 
ments  have  been  made  of  shock  waves  generated  by  large  pentolite  spheres  under 
river  or  brackish  water  at  distances  greater  than  12  charge  radii.  R.  S.  Price 
of  NSWC  summarizes  the  result  with  the  expression 

P  =  2.46  x  104  (W1/3/R)1*194  ,  (1) 

or 


where 

P  =  peak  pressure  in  pounds  per  square  inch.  Equation  (1) 

P  =  peak  pressure  in  kilobars.  Equation  (la) 

W  =  charge  weight  in  pounds 
R  =  radial  distance  in  feet 
Rq  =  charge  radius. 

3 

The  conversion  of  Equation  (1)  to  (la)  is  made  using  a  density  of  1.63  g/cm 
for  pentolite.  Corrections  for  salinity  and  temperature  variations  are  being 
made  to  the  data  on  which  Equation  (1)  is  based  so  that  an  equation  for  shocks 
in  fresh  water  for  the  large-scale  tests  will  be  available  in  the  future. 

The  purpose  of  this  report  is  to  make  available  to  others  the  data  we  have 
gathered  from  small-scale  tests  in  fresh  water  and  the  equations  that  were 
developed  which  closely  fit  the  data.  The  scalability  of  peak  pressure  and 
pulse  widths  as  functions  of  relative  radius,  R/R0,  is  a  significant  result. 

A  useful  analytical  expression  is  presented  which  accurately  reproduces  the 
pressure-vs-R/R0  data  over  the  range  of  R/R0  from  1.2  to  12.6.  In  the  limit, 
R/R0  =  1,  it  gives  a  pressure  which  agrees  with  the  value  obtained  by  Coleburn6 
using  the  Aquarium-Test  technique. 7  In  the  limit  R/R0  =  ~,  it  gives  zero 

6Coleburn,  N.  L.,  Chapman-Oouget  Pressures  of  Several  Pure  and  Mixed  Explosives, 
NOLTR  64-58,  Jun  1964. 

7Holton,  W.  C.,  The  Detonation  Pressure  in  Explosives  as  Measured  by  Transmitted 
Shocks  into  Water,  NAVORD  Report  3968,  Dec  1954. 


P  = 


18.88 


(R/R0) 


1.194 


(la) 
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pressure.  The  shock  pressures  and  time  histories  of  shocks  in  water  from  pentu 
lite  spheres,  calculated  by  Sternberg  and  Walker8,  are  in  fair  agreement  to  the 
results  of  this  study.  Also  there  is  reasonable  agreement  between  their 
calculations,  those  of  Rogers9,  and  the  far-field  results  given  by  Equation  (1) 


^Sternberg,  H.  M.  and  Walker,  W.  A.,  "Calculated  Flow  and  Energy  Distribution 
Following  Underwater  Detonation  of  a  Pentolite  Sphere,"  The  Physics  of  Fluids 
Vol .  14,  No.  9,  Sep  1971,  pp.  1869-1878. 

^Rogers,  P.  H.',  "Weak-Shock  Solution  for  Underwater  Explosive  Shock  Waves," 
Journal  of  the  Acoustical  Society  of  America,  Vol.  62,  1977,  pp.  1412-1419. 
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EXPERIMENTAL  PROCEDURE 


The  shock  velocity-versus-di stance  data  contained  in  this  report  is  derived 
mostly  from  Underwater  Sensitivity  Tests*,  the  experimental  arrangement  being 
shown  in  Figure  1.  The  spherical  charge  assembly  is  shown  in  Figure  2.  The' 
donor  is  an  82.2-mm  diameter,  cast  pentolite  sphere  (50%  TNT,  50%  PETN),  weighing 
472  ±  2  g;  density  =  1.63  g/cm3.  The  detonator  is  of  the  exploding  bridgewire 
type,  an  RP-80,  manufactured  by  Reynolds  Industries,  Inc.  It  is  7.11  mm  in 
diameter  and  fits  into  a  46-rnm  deep  hole  cast  into  the  sphere.  A  7.0 -mm  diameter 
by  9.5  mm-long  pellet  (density  =  1.6  g/cm3)  of  peutolite  is  placed  in  the  hole 
ahead  of  the  detonator  to  ensure  a  rapid  build-up  to  detonation  in  the  main 
charge.  The  space  around  the  detonator  leads  within  the  hole  is  filled  with  C-4 
plastic  explosive.  A  sealant  is  used  to  cover  the  connection  of  the  plastic 
sheath,  containing  the  detonator  leads,  to  the  pentolite  sphere. 

The  shocks  in  water  are  observed  by  a  Jacobs  high-speed  framing  camera, 
described  in  Appendix  A.  The  plane  of  focus  is  at  right  angles  to  the  axis  of 
the  detonator.  Diffuse,  reflected  back-lighting  is  provided  by  an  argon  explosive 
flashlamp  which  illuminates  a  white  cardboard  in  back  of  the  tank.  (An  exploding 
wire  light  source  with  a  Fresnel  lens  also  has  been  used  to  back-light  the  tank. 
However,  the  system  is  overly  sensitive  to  minor  disturbances  in  the  water  which 
tend  to  obscure  the  silouetted  outlines  of  the  acceptors.) 

A  calibration  of  the  underwater  system  is  obtained  from  the  framing-camera 
records.  For  reasons  of  economy  and  expendiency,  the  Underwater  Sensitivity  Test 
experiments  and  a  calibration  of  the  system  are  carried  out  simultaneously. 

A  typical  sequence  of  frames  is  shown  in  Figure  3.  The  interframe  time  is 
3.63  us.  (Two  framing  camera  film  strips  are  exposed  in  each  shot,  the  frames 
of  the  two  tracks  being  separated  by  1.81  ys.)  Not  all  of  the  records,  or 
portions  of  records,  are  used  in  determining  the  calibration.  On  the  records 
from  some  experiments,  one  or  more  acceptors  are  fairly  close  to  the  donor, 
thus  limiting  the  number  of  x-t  readings  that  are  made  in  their  direction.  On 
other  records,  the  degree  of  blur  due  to  motion,  or  to  being  somewhat  out  of 
focus,  is  enough  to  eliminate  al'  ,  or  part,  of  the  record  from  the  calibration 
analysis.  For  the  reasons  given  above,  acceptable  readings  usually  are  obtained 
in  only  one  or  two  of  the  four  major  directions  on  the  record,  i.e.,  toward  the 
top,  bottom,  left,  or  right.  Of  these  four  directions,  the  first  two  usually 
produce  the  sharpest  images  on  the  film.  As  a  result,  only  one  record  is 
included  in  the  calibration  with  the  snock-front  image  moving  at  right  angles 
to  the  direction  of  slit  motion. 

♦Sometimes  called  the  Aquarium  Test,  but  should  not  be  confused  with  the  underwater 
test  of  reference  7  from  which  detonation  pressures  are  obtained. 
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The  addition  to  photographic  coverage*  lithium  niobate  gages  are  used  in 
some  experiments  to  measure  peak  pressures  and  time  histories  at  several  distances 
from  the  pentolite  sphere.  The  lithium  niobate  rage  has  the  advantage  of  being 
made  from  commercially  grown  synthetic  crystals  and  its  calibration  apparently 
is  independent  of  manufacturer. 10 

The  lithium  Riobate  gage  is  used  with  capacitive  loading  as  shown  in 
Figure  4.  The  charge  on  the  capacitor,  generated  by  hydrostatic  pressure,  can 
be  recorded  by  an  oscilloscope.  The  lithium  niobate  hydrostatic  piezoelectric 
constant  of  6.3  x.10"®  couiomb/cm^-kbar  is  used  to  convert  charge  to  pressure. 

The  oscilloscope  input  impedance  must  be  large  so  that  the  discharge  time  of 
the  capacitor  through  the  scope  impedance  is  small  compared  to  the  time  of  the 
event.  Another  consideration  is  to  load  the  gage  with  a  capacitance  that  is 
nearly  100  times  that  of  the  gage  cable  so  that  any  changes  in  cable  capacitance 
due  to  shock  loading  can  be  ignored. 

Typical  gage  records  are  shown  in  Figure  5.  There  exists  an  extra  oscilla¬ 
tion  superimposed  on  the  voltage  signal  due  to  the  pressure  pulse.  This  extra 
signal  is  due  to  a  normal  mode  oscillation  of  the  lithium  niobate  crystal.-1- 
In  reading  the  records,  these  normal  mode  oscillations  are  accounted  for  by 
drawing  a  line  symmetrically  through  them.  The  peak  pressure  is  determined 
by  the  intersection  of  the  line  through  the  normal  mode  oscillations  and  a 
projected  line  for  the  rise  time  of  the  pressure  pulse. 


^Graham.  R.  A.,  "Pressure  Dependence  of  the  Piezoelectric  Polarization  of 
LiNbOg  and  LiTa03,"  Ferroelectrics,  Vol .  10,  1976,  pp.  65-69. 

^Goldstein,  S.,  and  Johnson,  J.  N. ,  "Aquarium  Tests  on  Aluminized  ANFO," 
Pre-prints  of  the  Seventh  Symposium  on  Detonation,  16-19  Jun  1981,  pp.  524-530. 
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TREATMENT  OF  DATA 


The  pressure  in  the  water,  P,  as  a  function  of  distance,  x  or  R/R0,  is 
derived  front  experimentally  determined  values  of  the  shock  velocity  as  well  as 
from  pressure  gages.  The  derivation  from  the  shock  velocity  is  made  using  the 
momentum  relation,  P  «  pUsu,  and  the  Rice-Walsh  equation  of  state  for  water, 

Us  *  c0  +  10.990  ln(l  ♦  u/5i„9}.i2  Combining  the  two  equations  results  in 


P  » 


51.9  pUs 


107990 


where 


P  *  pressure  (kba.') 

Us  =  shock  velocity  (mm/ys) 
u  *  particle  velocity  (mm/ys) 
cQ  a  sound  velocity  (1.483  mm/ys) 
p  =  density  (0.998  g/cm3). 


To  obtain  shock  velocities  from  the  framing  camera  records,  the  x-t  data 
are  plotted,  a  smooth  curve  is  drawn  through  the  data,  and  the  slopes  are  drawn 
at  given  values  of  x  from  10  to  120  mm.  However,  before  the  x-t  data  ere  plotted 
certain  corrections  are  made  both  in  time  and  in  distance  (position).  The 
position  correction  is  a  consequence  of  shock-front  motion  during  the  exposure 
time  of  the  camera  slits  (focal  plane  shutters).  If  the  velocity  is  constant, 
the  position  error  is  the  same  for  all  frames.  However,  any  difference  in 
shock  velocity  occurring  between  frames  results  in  position  errors  that  affect 
velocity  calculations.  For  the  detonated  sphere  the  shock  velocity  initially 
is  high,  but  rapidly  decaying.  This  results  in  fairly  large  errors  in  position 
very  close  to  the  charge.  The  errors  in  position  become  relatively  small  a 
short  distance  away  from  the  charge  surface.  A  further  discussion  of  correction 
due  to  motion  blur  is  in  Appendix  A. 

As  to  the  time  correction,  the  direction  of  shock  motion  on  the  film, 
relative  to  the  camera  slits,  changes  the  interframe  time  if  it  is  in  the  same 

12Rice,  M.  H.  and  Walsh,  J.  M. ,  "Equation  of  State  of  Water  to  250  Kilobars," 

The  Journal  of  Chemical  Physics,  Vol .  26,  No.  4,  Apr  1957,  pp.  824-830. 
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direction  as  the  slit  motion,  or  if  it  opposes  it.  The  time  corrections  are 
added  if  the  slits  and  image  are  moving  in  the  same  direction  and  subtracted 
if  they  oppose  each  other.  No  time  correction  is  necessary  if  the  shock-front 
image  motion  is  at  right  angles  to  the  motion  of  the  slits.  This  also  is 
discussed  in  Appendix  A.  The  uncorrected  x-t  data  are  found  in  Table  A-l  and 
the  corrected  data  in  Table  A-3  of  Appendix  A. 


The  shock-velocity  data,  Table  1,  are  derived  from  eight  Underwater 
Sensitivity  Test  experiments,  having  14  usable  x-t  framing  camera  sequences* 
over  the  range  from  x  *  10  to  60  mm.  Data  for  seven  shots  (13  sequences)  are 
obtained  from  x  -  10  to  70  or  80  nm,  five  shots  (11  sequences)  from  x  =  10  to 
100  mm,  and  three  shots  (5  sequences)  from  x  =  10  to  120  mm.  In  Table  1,  two 
rows  of  instantaneous  velocity  values,  Us,  are  listed  for  each  value  of  x. 

The  upper  row  gives  velocities  obtained  from  all  of  the  x-t  sequences  up  to 
fourteen.  The  lower  row  gives  the  results  from  the  eight  individual  experiments 
(shots).  Where  ware  than  one  x-t  sequence  is  used  per  shot,  the  velocity  values 
are  averaged.  The  experimental  values  of  the  shock  velocity  (eight  shots),  Us, 
the  average  deviation  of  the  mean,  a,  and  the  number  of  measurements  made  are 
listed  in  Table  2  for  various  values  of  x  and  R/R0.  Two  shock  velocities 
obtained  from  smear-camera  records  also  are  included  in  Table  2.  The  peak 
pressures  in  the  Table  are  obtained  from  Equation  (2)  and  the  experimental 
values  listed  for  Us« 


The  best  fit  to  the  data  over  the  entire  range  is  obtained  using  the 
following  equation,  the  coefficients  being  determined  by  a  least  squares 
technique, 


Us  "  co 


(3) 


where 


Rx  =  R/Rq 
A  *  -0.406213 


B  =  1.30135  . 

This  equation  gives  a  limit  of  Us  *  6.053  mm/ys  at  Rj  =  1  which  is  adjusted  to 
agree  with  the  experimental  value  obtained  by  Coleburn®  for  50/50  pentolite. 
Also,  It  converges  to  sonic  velocity  in  the  limit,  Ri  s  ®.  Several  different 
expressions  were  tried  before  settling  on  Equation  (3).  The  varlow*  cxprassions 
attempted  and  their  fits  to  the  data  are  presented  in  Appendix  B. 

*In  eight  shots  with  16  film  strips,  the  number  of  possible  x-t  sequences  (up, 
down,  left,  and  right)  is  64. 
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TABLE  1.  INSTANTANEOUS  SHOCK  VELOCITIES  DERIVED  FROM 
FOURTEEN  FRAMING- CAMERA  SEQUENCES  (EIGHT  SHOTS) 


•j  \ 
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TABLE  2.  SUMMARY  OF  THE  MEASURED  INSTANTANEOUS  SHOCK  VELOCITIES  AND 
THE  CORRESPONDING  PRESSURES  OBTAINED  FROM  EQUATION  (2) 


X 

(mm) 

R/R0 

NUMBER 

MEASURED 

(mm? us) 

AV.  DEV. 

OF  MEAN 

P  a 
(kbar)3 

10 

1.243 

8 

4.405 

0.127 

69.5 

20 

1.487 

8 

3.649 

0.070 

41.2 

30 

1.730 

8 

3.157 

0.038 

26.9 

40 

1.973 

8 

2.829 

0.028 

19.1 

50 

2.217 

8 

2.596 

0.014 

14.3 

60 

2.460 

8 

2.434 

0.019 

11.4 

70 

2.703 

7 

2.311 

0.027 

9.4 

80 

2.946 

7 

2.218 

0.019 

8.0 

100 

3.433 

5 

2.084 

0.015 

6.1 

120 

3.920 

3 

1.975 

0.005 

4.7 

89 

8.02b 

1 

1.686c 

0.03 

•  1.63 

143 

12 . 27b 

1 

1 . 583c 

0.02 

0.75 

aFrom  Equation  (2) 

bRQ  c  12.7  mm,  all  others  are  41.1  mm 

c$mear-camera  measurements  of  shock  velocity;  the  other  velocities  listed  are 
from  Jacobs  framing-camera  records.  Gage  measurements  also  were  made  in  these 
two  experiments. 
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The  peak  pressures  and  pulse  half-widths  are  measured  by  lithium  niobate 
gages  for  a  number  of  cast  and  pressed  pentolite  spheres  with  radii  of  L27, 
4.11,  and  6.84  cm.  Information  on  the  spheres  is  given  in  Table  3.  The  peak 
pressure  and  pulse  widths  at  one-half  peak  pressure  are  given  in  Table  4. 

The  measured  pressures  span  a  range  of  0.8  to  5.9  kbar  for  a  relative  distance, 
R/Rp,  from  12.3  to  3.5,  respectively.  The  peak  pressures  for  the  same  R/R0  in 
Table  4  are  the  same  within  experimental  error.  The  pulse  half-widths  ranged 
from  10  to  66  ys.  The  pulse  half-widths  divided  by  R0  give  the  same  results 
for  the  different  size  spheres  when  the  measurements  are  made  at  the  same  R/R0. 
The  one  exception  in  Table  4  is  for  t/R$  *  11.73  ps/cm  at  R/Rq  *  8.01.  No 
adequate  explanation  exists  for  its  departure  from  the  noted  trend. 
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TABLE  3.  PENTOLITE  SPHERE  AND  LITHIUM  NIOSATE  GAGE  SPECIFICATIONS 


SHOT 

NO. 

GAGE 

NO. 

CHARGE  RAD. 
(cm) 

DENSITY 

(g/cm3) 

GAGE  CAP.3 
(10-9F) 

TOTAL  CAP. 
(10-9F) 

GAGE  DI$T.C 
(cm) 

1 

1 

1.27 

1.61b 

0.45 

10.3 

15.56 

2 

1 

1.27 

1.62b 

0.45 

10.3 

10.16 

3 

1 

4.11 

1.65 

— 

40.0 

18.06 

4 

1 

4.11 

1.64 

— 

79.0 

14.26 

4 

2 

4.11 

1.64 

-- 

40.0 

20.41 

5 

1 

4.11 

1.62 

1.5 

121.0 

32.94 

5 

2 

4.11 

1.62 

1.5 

81.0 

50.47 

6 

1 

6.84 

1.62 

"2.0 

121.0 

44.49 

6 

2 

6.84 

1.62 

■  "2.0 

81.0 

54.82 

*5 

Includes  capacitance  of  cable 
^Pressed 

cFrom  center  of  sphere 
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TABLE  4.  GAGE  MEASUREMENTS  OF  SHOCKS  FROM 
THREE  SIZES  OF  PENTOLITE  SPHERES 


SHOT 

NO. 

R/R0 

PEAK 

PRESSURE 

(kbar) 

t, PULSE 
HALF-WIDTH 
(us) 

t/R0 

(ws/cm) 

WATER 

TEMP. 

(°C) 

1 

12.27 

0.803  ±  0.030 

15.23  ±  0.61 

12.01  ±  0.48 

21.1 

2 

8.02 

1.447  ±  0.043 

12.46  ±  0.93 

9.83  ±  0.58 

22.2 

3 

4.39 

3.823  ±  0.126 

24.56  ±  0.85 

5.98  ±  0.21 

17.8 

4 

3.47 

5.862  ±  0.251 

10.32  ±1.93 

2.51  ±  0.47 

17.8 

4 

4.97 

3.048  ±  0.080 

24.66  ±  1.16 

6.00  ±  0.28 

17.8 

5 

8.01 

1.543  ±  0.076 

48.19  ±  2.96 

11.73  ±  0.72 

20.6 

5 

12.28 

0.817  ±  0.033 

50.08  ±  2.42 

12.18  ±  0.59 

20.6 

6 

6.50 

2.038  ±  0.075 

57.15  ±  2.87 

8.36  ±  0.42 

26.0 

6 

8.01 

1.449  ±  0.061 

65.88  ±  1.48 

9.63  ±  0.22 

26.0 
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DISCUSSION 


A  useful  form  of  the  data  is  in  terms  of  peak  pressure  as  a  function  of 
distance  in  water  from  the  pentolite  sphere.  Combining  Equation  (3)  with 
Equation  (2)  results  in 


P  =  76.814 


1  + 


1.8298 


1.3014 


'1 


-0.4062 


exp 


0.2469 


R,1-3014-!). 


5062 


where 


P  =  peak  pressure  in  kilobars 
R1  =  relative  distance,  R/RQ,  from  1,24  to  14.3. 


(4) 


The  coefficients  given  here  differ  somewhat  from  those  given  in  reference  13.  • 
The  change  is  due  to  the  inclusion  of  more  refined  x-t  data  and  velocity 
measurements  in  the  later  calculations.  In  Figure  6,  the  solid  line  gives 
the  peak  pressure  as  a  function  of  R/R0  as  obtained  from  Equation  (4).  The 
lithium  niobate  data,  shown  as  bars,  are  in  good  agreement  with  the  solid 
line.  The  dashed  line  in  Figure  6,  representing  Equation  (1)  is  in  fair 
agreement  with  the  data  for  pressures  below  five  kilobars.  Note  that  Equation 
(4)  predicts  pressures  slightly  lower  than  Equation  (1)  for  large  values  of 
Rl,  i.e.,  R/Ro>10.  For  the  most  part,  the  small  differences  in  the  predicted 
results  from  these  two  equations  for  large  R/R0  are  likely  to  be  due  to  gage 
inaccuracies  and  to  Equation  (1)  being  based  on  data  from  large-scale  field 
tests  in  brackish  water  with  variations  in  temperature.  The  pressure  calculated 
from  the  measured  shock  velocity  at  R/R0  =  8.02  in  Table  2  is  11%  higher  than 
the  gage  measurement,  listed  in  Table  4  for  the  same  experiment.  A  velocity  of 
only  1%  lower  than  shown,  however,  results  in  a  pressure  that  agrees  with  the 
corresponding  gage  measurement  to  within  experimental  error. 


The  measured  pulse  half-widths  as  a  function  of  R/R0  are  given  in  Figure  7. 
The  gage  pulse  width  data  are  in  fair  agreement  with  the  calculated  results,  the 
measured  pulse  half-widths  being  higher  than  the  code  calculation  for  R/R0  values 
greater  than  seven.  For  R/R0  less  than  five,  the  measured  half-widths  are  less 
than  calculated.  The  record  at  3.5  has  a  structure  to  it  which  other  records 
did  not  show.  It  is  possible  that  this  gage  w  s  breaking  under  high  loading. 

This  possibility  is  further  supported  by  the  f«Jt  that  a  gage  at  R/R0  3  1.89 
gave  no  record  (not  included  in  Table  4)  and  we  suspect  that  its  elastic  yield 
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FIGURE  6.  PEAK  PRESSURE  VERSUS  RELATIVE  DISTANCE  FROM  PENTOLITE  SPHERES 
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point  was  exceeded  causing  irreversible  damage  to  the  crystal.  The  trend  in 
the  data  is  in  agreement  with  the  hydrodynamic  code  calculations. 8 

Shock  data  for  distances  greater  than  R/R0  =  12  are  of  secondary  interest 
to  this  work.  However,  to  make  this  a  more  complete  report,  the  work  of  Sternberg 
and  Walker^  and  that  of  Rogers^  will  be  briefly  discussed.  The  hydrodynamic 
code  calculations  by  Sternberg  and  Walker^  are  in  good  agreement  with  the  field- 
test  data  represented  by  Equation  (1).  The  paper  by  Rogers*  successfully  derives 
an  analytical  expression  for  peak  pressure  and  pulse  width  for  exponentially 
decaying  shocks  as  a  function  of  distance  from  the  charge.  To  use  these  expres¬ 
sions  requires  knowing  the  shock-pulse  amplitude  and  width  at  one  position. 

The  theory  then  predicts  the  pulse  amplitude  and  width  at  later  times  and 
positions.  This  theory  was  not  expected  to  give  accurate  results  for  shock 
pressures  of  greater  than  about  two  kilobers.  In  Appendix  C  are  shown  the 
results  of  fitting  the  peak  pressures  and  time-pulse  widths  obtained  in  this 
study  to  the  analytical  expressions  derived  by  Rogers. 
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CONCLUSIONS 


Shock  wave  velocities,,  peak  pressures  and  pulse  widths  at  one-half  peak 
pressure  have  been  measured  as  a  function  of  distance  in  fresh  water  from 
pentolite  spheres.  An  analytic  function  that  reproduces  the  pressure  data 
within  experimental  error  over  the  relative  distance  (R/Ro)  interval  from  1.2 
to  12  has  been  found.  This  expression  is  empirical  and  results  from  the 
observation  that  the  shock  velocity  minus  the  sound  velocity  is  proportional 
to  X/fRB  +  A).  The  pressure  equation  goes  to  the  desired  limits  of  162  kbar 
at  R/R0  =  1  (the  experimental  value  determined  by  Coleburn)  and  0  kbar  at 
R/R0  =  «.  The  measured  peak  pressures  agree  within  experimental  error  with 
Equation  (4). 

The  calculated  shock  pressures  and  time  histories  of  shocks  in  water  from 
pentolite  spheres  by  Walker  and  Sternberg  are  in  fair  agreement  with  the  present 
experimental  results.  Extrapolation  of  the  present  data  by  Equation  (4)  to 
R/R0  >  13  results  in  pressures  slightly  lower  than  measure  in  large-scale  field 
tests. 

The  pressure  pulse  in  water  from  the  pentolite  sphere  can  be  scaled  as  a 
function  of  the  sphere  radius.  This  indicates  that  most  or  all  of  the  pentolite 
detonates  even  for  very  small  spheres  of  1.27  cm  radius.  This  is  consistent 
with  the  known  short  run  distance  to  detonation  for  pentolite. 
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APPENDIX  A 

CORRECTION  OF  DISTANCE-TIME  DATA  DERIVED 
FROM  JACOBS  CAfCRA  RECORDS 


THE  JACOBS  CAMERA 


The  Jacobs  .camera  records  a  total  of  six  rows  of  24  x  18  mm  pictures,  three 
rows  on  each  of  two  70  m  film  tracks.  A  maximum  of  216  pictures  are  available. 

In  a  focal  plane  shutter  camera  only  a  narrow  band  of  the  frame  is  exposed  by 
the  moving  aperture  (a  slit  in  this  case)  at  a  single  instant  of  time.  Ultimately 
the  entire  field  of  view  is  scanned  as  the  slit  moves  across  the  whole  frame. 

In  this  camera  six  pictures  are  partially  exposed  simultaneously  (one  per  row) 
through  six  slits,  the  slits  being  displaced  one-sixth  of  a  frame  apart.  Thus, 
for  the  duration  of  the  recording  portion  of  the  camera's  duty  cycle,  the  field 
is  continuously  scanned  by  six  narrow  bands  over  and  over  again,  each  scan  by 
each  band  exposing  a  separate  picture.  The  shutter  action  is  illustrated  in 
Figure  A-l. 

In  this  camera  time  on  the  film  track  increases  in  two  modes,  between  frames 
and  within  frames”.  Within  a  single  frame  the  focal  plane  shutter  (actually  the 
1  :  1  image  of  .the  slit)  sweeps  across  the  frame  area'exposing  the  film  one  slit 
width  at  a  time.  The  slit  widths  can  be  set  in  either  track  at  0.002,  0.005, 
0.010,  0.020,  0.040,  and  0.080  inch.  At  a  rotating  mirror  speed  of  300  rps  the 
focal  plane  shutter  sweep  velocity  is  2.298  bsr/vs.  (This  is  the  speed  used 
throughout  the  Underwater  Sensitivity  Test  experiments.)  Thus,  time  within  a 
frame  increases  at  a  rate  of  0.435  ys/mm.  At  this  rate  the  slit  traverses  the 
25  mm  distance,  from  a  point  in  one  frame  to  the  corresponding  point  in  the 
next  frame,  in  10.88  ys.  This  is  the  interframe  time  for  horizontally  adjacent 
images.  Since  these  shutters  move  simultaneously  along  each  track,  phased  one- 
sixth  of  a  frame  apart,  time  advances  by  3.628  ys  between  frames  in  all  three 
rows  in  each  track  and  1.814  ys,  alternating  between  frames  in  both  tracks. 

For  further  details  concerning  the  Jacobs  camera,  see  reference  A-l. 

CORRECTIONS  FOR  MOTION  BLUR 

In  reading  the  framing  camera  records,  the  shock-front  image  position  on  the 
film,  including  the  motion  blur,  ax',  is  denoted  by  x'  +  ax'.  The  position 
correction,  however,  is  made  after  the  original  readings  are  converted  to  full 
scale,  i.e.,  they  are  divided  by  the  magnification  factor  (actually  a  reduction 
factor),  M,  giving  values  of  x  +  ax.  The  real  (corrected)  shock-front  position, 
x,  is  obtained  by  subtracting  out  ax,  after  ax  has  been  determined.  The 


^Jacobs,  S.  J.,  McLanahan,  J.  0.,  Jr.,  and  Whitman,  E.  C.,  "A  High-speed 
Focal -Plane  Shutter  Framing  Camera,"  Journal  of  the  Society  QfJlQMpn 
Picture  and  Television  Engineers,  Vol .  72  No.  12,  Dec  1963,  pp.  923-926. 
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RAMING  SEQUENCE 
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uncorrected  readings  of  x'  +  ax'  from  the  records  are  given  In  Table  A-l  for 
fourteen  sequences,  along  with  the  corresponding  (full  scale)  values  of  x  ♦  ax. 

In  most  cases,  M  can  be  considered  constant  throughout  the  sequence.  However, 

In  a  few  cases  M  changes  gradually  over  the  distance  being  read  over  the  record. 
For  the  most  part,  this  is  due  to  some  tilt  of  a  mirror  placed  in  the  optical 
path  within  the  firing  chamber  to  give  a  greater  distance  to  the  object.  The 
image  of  a  graduated  scale  on  the  film.  In  the  plane  of  focus  within  the  tank, 
permits  a  check  of  the  constancy  of  the  magnification  factor.  The  values  of  M 
used  are  not  listed  in  the  Table,  but  can  be  determined  by  dividing  x  +  ax  by 
x'  +  ax'. 

How  a  correction  is  made  for  image  motion  blur  depends  on  the  direction 
taken  by  the  point  being  followed  on  the  shock-front  image.  For  a  slit  width, 

W,  and  a  slit  velocity,  Ur,  the  correction  for  image  blur,  ax',  Is  different 
for  each  of  the  three  following  conditions:  Case  1,  when  the  shock-front  image 
motion  Is  in  the  same  direction  as  the  slit  motion,  i.e.,  downward;  Case  2, 
when  the  shock-front  motion  is  in  the  opposite,  upward  direction;  and  Case  3, 
when  the  shock-front  motion  is  at  right  angles  to  the  direction  of  slit  motion. 

In  Case  1,  shock  and  slit  in  the  same  direction,  the  effective  slit  width 
Is  W  +  ax',  ax'  being  the  amount  of  blur  on  the  film  and  is  incorporated  in  the 
distance  value  read  from  the  framing  camera9  i.e.,  x'  +  ax'.  This  is  illustrated 
on  the  right  side  of  Figure  A-2.  At  t  a  ti  the  shock-front  Image  just  overtakes 
the  leading  edge  of  the  slit  image,  Aj.  The  exposure  of  the  film  lasts  until  the 
back  of  the  slit,  Bi,  overtakes  the  shock-front  image  at  t  *  t2  (position  B2  in 
the  sketch.)  At  the  Instant  the  shock-front  image  is  cut  off,  the  leading  edge 
of  the  slit  has  moved  from  Ai  to  A2  or  a  distance  equal  to  the  slit  width,  W, 
plus  the  motion  blur,  ax'.  Thus,  the  exposure  time  of  the  moving  shock-front 
image,  z,  is  (W  +  Ax*)/Uf.  In  Case  2,  with  shock  and  slit  moving  in  opposing 
directions,  the  effective  slit  width  is  W  -  ax'.  This  can  be  deduced  from  the 
slit  and  shock-front  positions  shown  at  the  left  side  of  Figure  A-2.  (Note  that, 
in  reality,  the  image  of  the  expanding  shock  front  cannot  be  a  true  circle  and 
tj  cannot  have  the  same  value  for  all  three  cases  on  the  record.  The  shock 
fronts  are  drawn  spherically  here  for  simplicity  and  this  in  no  way  affects 
the  results  when  each  case  Is  considered  separately.)  It  can  be  shown  that  z' 
also  is  equal  to  ax'/(W‘D‘s),  0$  being  the  average  velocity  of  the  shock-front  _ 
image  across  the  effective  slit  width.  Therefore,  z'  *  (W  ±  Ax')/llf  *  ax'/(M*Us). 
Solving  for  blur  on  the  film, 


W-M-U 

AX'  *  u^rwv;  ’ 

where  the  sign  in  the  denominator  is  minus  when  Uf  and  M*US  move  in  the  same 
direction  and  positive  when  they  move  in  opposition.  It  follows  then  that  the 
full-scale  (object)  motion  blur  is 


AX 


AX1 

M 


(A-l) 


L£  A-l.  UNCORRECTED  x-t  DATA  FROM  EIGHT  EXPERIMENTS  HAVING 
FOURTEEN  USABLE  SEQUENCES  (x'  Ax'  and  x  +  Ax  In  ram) 


m 


-  .U  -  M 


In  Case  3,  with  the  shock-front  image  motion  at  right  angles  t:  the  direction 
of  slit  motion,  the  exposure  time  is  taken  to  be  equal  to  the  slit  duration.  This 
is  indicated  in  the  center  of  Figure  A-2.  The  slit  duration  is  defined  here  as 
the  time  it  takes  the  slit  to  travel  its  own  width,  i.e.,  W/Uf.  The  motion  blur 
on  the  film,  ax',  is  defined  as  (M- W)/Us  in  this  case,  and  the  object  motion  blur, 
ax,  is  W/0$. 

In  determining  the  correction  curve,  a  series  of  values  of  x  (the  true,  full- 
scale,  shock-front  position}  are  chosen  at  convenient  intervals  from  0  to  120  mm. 
The  required  values  of  Oc  are  determined  in  the  following  way.  First,  the  instan¬ 
taneous  shock  velocity,  us,  is  found  for  each  value  of  x  selected,  using 
Equation  (3).  The  values  of  corresponding  to  the  true  position,  x,  are  listed 
in  the  first  two  columns  of  Table  A-2.  An  approximate  value  of  ax  for  Cases  1  or 
2  is  obtained  by  substituting  Us  for  CL  in  Equation  (A-l).  The  value  of  M  used 
here  is  25.3,  the  mean  yalue  of  all  x-t  sequences.  (Actually,  one  could  use  any 
value  of  M  from  24  to  26  without  affecting  the  results  appreciably.)  In  Case  3, 
the  substitution  is  made  in  the  expression  ax  =  w/0$.  The  approximate  values  of 
ax,  so  determined,  are  added  to  x  and  U$  is  calculated  at  that  position.  The  mean 
between  the  two  values  at  x  and  x  +  ax  is  used  to  determine  a  more  refined  value 
of  ax.  The  process  is  repeated  to  obtain  the  final  value  of  0S  and  ax  given  in 
Table  A-2.  The  resulting  correction  curves  for  five  combinations  of  two  slit 
widths  and  three  directions  of  slit  motion  are  shown  in  Figure  A-3.  Using  correc¬ 
tions  taken  from  the  curves,  the  uncorrected  data  in  Table  A-l,  x'  +  ax'  for  the 
image  and  x  +  ax  ^or  full  scale,  are  converted  to  the  corrected  distance,  x,  given 
in  Table  A-3.  The  corrected  time,  t,  also  is  Included  and  is  discussed  next. 

CORRECTIONS  FOR  TIME 

After  the  values  of  x',  i.e.,  x/M,  are  determined,  the  correction  for  time, 
At,  can  be  made  simply  by  dividing  a'  by  the  slit  velocity,  U|r  (2.298  mm/ys),  i.e. 
At  =  x'/Uf.  If  the  shock  front  and  slits  are  moving  in  the  same  direction.  At  is 
added  to  the  accumulated  time,  t'.  In  each  track  the  uncorrected  time,  t' , 

increases  by  10.883  ys  per  frame  along  a  row  and  by  3.628  ys  when  taking  into 

account  both  rows  and  columns.  If  the  direction  of  shock  motion  is  in  opposition 

to  the  slit.  At  is  subtracted  from  t' .  There  is  no  correction  when  the  shock 

motion  is  at  right  angles  to  the,  slit  motion.  The  derivation  of  the  time  correc¬ 
tion,  At,  can  be  demonstrated  by  Figure  A-4.  Three  consecutive  frames,  running 
parallel  to  the  film  edge,  are  shown  (upper  row)  in  which  an  element  of  the  shock- 
front  image  is  advancing  in  the  same  direction  as  the  focal  plane  shutter.  The 
position  of  the  shock-front  element  in  the  first  frame  is  taken  as  zero  reference, 
(t0,  Xg).  The  primes  added  to  t«  in  succeeding  frames  denote  zero  time  plus 
accumulated  interframe  times.  The  resulting  reference  time  lines,  t0,  t q  and  t^' , 
are  in  the  same  relative  position,  x0,  within  each  frame.  It  is  seen  that  in 
order  for  the  shock  front  to  be  exposed  in  the  next  frame,  the  slit  must  move  from 
a  position,  x0,  to  xi  which  lies  to  the  right  of  the  reference  time  line,  to  Thi 
means  that,  in  effect,  the  time  between  frames  increases  as  far  as  the  shock-front 
image  is  concerned.  The  correction  for  time,  t,  is  the  accumulated  increase  in 
distance  on  the  film,  xi  -  x0,  x3  -  x0,  x3  -  x0»  etc.,  divided  by  the  slit 
velocity,  Uf.  For  shock-front  image  motion  in  the  opposite  direction  (lower  row 
of  frames),  it  is  seen  that  the  time  between  frames  decreases  for  the  shock-front 
image. 
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TABLE  A-2.  THE  INSTANTANEOUS  SHOCK  VELOCITY,  Us,  THE  AVERAGE  SHOCK 
VELOCITY,  Us,  AND  MOTION  CORRECTION,  Ax,  AS  A  FUNCTION  OF  x 
FOR  TWO  SLIT  WIDTHS  AND  THREE  SHOCK  DIRECTIONS 
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-3.  CORRECTED  x-t  DATA  FROM  EIGHT  EXPERIMENTS  HAVING 
FOURTEEN  USABLE  SEQUENCES  (  x  in  mm,  t  in  ys) 
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In  the  preceding  paragraph  it  is  assumed  that  the  spacing  between  the  position 
of  x0  on  adjacent  frames  is  always  the  same.  In  reality  this  is  not  so  since  minor 
variations  in  frame  spacing  are  inherent  in  the  camera.  In  the  analysis  in  this 
report  these  variations  make  no  significant  difference  in  fitting  a  curve  to  the 
x-t  data.  Therefore,  this  correction  is  ignored  in  our  analysis.  However,  the 
variation  in  frame  spacing  can  lead  to  relatively  large  time  errors  in  some 
applications,  especially  when  only  a  few  frames  are  available  in  a  sequence. 
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APPENDIX  B 

THE  SELECTION  OF  SEVERAL  ANALYTICAL  FORMS  FOR 
FITTING  SHOCK  VELOCITY-TO-DISTANCE  DATA 


i‘\- 


If  Us,  from  Table  2,  is  plotted  as  a  function  of  1/x,  Figure  B-l,  it  is  seen 
that  a  nearly  straight  line  can  be  drawn  through  most  of  the  datum  points  between 
x  *  40  and  120  mm.  The  equation  for  this  line  is 


U$  =  c’ 


1  + 


A* 


3-D 


where 


L:v 


c'  A'  -  the  slope  and 

c ’  =  the  Us-axis  intercept. 

The  values  of  c'  and  A'  are  determined  by  a  least  squares  fit  to  the  data  in 
Table  2.  In  the  range  of  x  *  50  to  100  mm,  the  value  of  c'  is  1.5762  nsn/ys  and 
A'  is  32.493  mm  with  a  quadratic  mean  error  (QME)  of  0.0025.  This  is  an  extremely 
good  fit.  From  x  *  40  to  120  mm,  c!  is  1.5739  and  A'  is  32.306  with  a  QME  of 
0.0153;  a  good  fit.  However,  for  x  less  than  40  mm  the  fit  is  not  good.  This  is 
expected,  of  course,  since  as  x  approaches  zero-  the  calculated  velocity  approaches 
infinity.  The  calculated  velocities,  using  Equation  (B-l),  are  compared  with  the 
measured  velocities  in  Table  B-l.  The  velocities  listed  correspond  to  data  in  the 
range  of  x  *  50  to  100  aw  used  to  calculate  c'  and  A'. 

Although  Equation  (B-l)  is  not  valid  for  x  much  less  than  40  mm,  it  can  be 
used  as  a  isaans  of  checking  all  of  the  x-t  data  sets  for  consistency  and  for 
adjusting  all  of  the  x-t  data  so  that  they  can  be  compared  in  a  single  x-t  plot. 
Equation  (B-l)  is  the  only  expression  out  of  the  several  considered  in  which  x 
can  be  calculated  with  reasonable  ease.  To  obtain  x,  the  Equation  is  first  put 
in  the  differential  form,  xdx/(x  +  A' )  =  c'dt.  On  integrating. 


K  *  c't  +  A'ln(x  +  A')  -  x. 


(B-2) 


Since  there  is  no  direct  algebraic  solution  for  x,  Newton's  Method  of  solution  to 
f (x)  -  0  is  used,  i .e. > 


f (x)  =  x  -  A'ln(x  +  A')  +  K  -  c't  »  0. 


(B-3) 


One  can  solve  for  the  constant  of  integration,  K,  in  Equation  (B-2)  by  using 
previously  determined  values  of  c'  and  A’  and  the  corrected  values  of  x  and  t  from 
Table  A-3,  Appendix  A.  It  is  found  that  K  is  fairly  constant  for  values  of  x 
ranging  from  30  to  125  mm.  The  mean  value  of  K,  i.e.,  R,  is  determined  over  this 
constant-value  range  for  each  x-t  set.  However,  the  value  of  R  is  observed  to 
vary  considerably  from  set  to  set.  This  occurs  because  the  values  of  t  are  taken 

B-l 
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to  be  evenly  spaced  interframe  times,  before  corrections  are  made,  and  x  =  0, 
t  =  0  can  fall  anywhere  between  the  frames.  Actually,  the  time  at  which  shock- 
front  motion  begins  can  not  be  pin-pointed  on  any  one  record.  At  the  start  the 
shock  velocity  is  about  6  mm/ys,  but  it  decays  rapidly.  Thus,  in  the  first  10  mm 
of  travel,  the  shock-front  image  is  recorded  in  only  one  frame,  or  none  at  all.  In 
fact,  out  of  the  fourteen  sequences,  there  are  only  six  shock-front  images  recorded  • 
within  the  first  10  nm  of  travel. 

All  of  the  x-t  data  can  be  put  in  proper  time  correspondence  by  subtracting 
the  value  of  t  found  at  x  =  0  in  Equation  (B-2)  from  the  values  of  t  listed  for 
each  x-t  set  in  Table  A-3,  Appendix  A.  When  the  adjustments  are  made,  the  value 
of  K  is  found  to  equal  113.1089  mm  for  all  of  the  x-t  sets.  The  adjusted  x-t  data 
are  plotted  in  Figure  8-2.  The  close-in  shape  of  the  x-t  curve  is  now  apparent 
and  the  actual  origin  of  the  x-t  curve  is  seen  to  be  located  at  about  t  =  -1.8  us. 

Finally,  the  corrected  values  of  t  from  Table  A-3  and  the  determined  values  of 
c1  and  A'  are  used  in  Equation  (B-3)  to  calculate  x.  The  results,  Xra]c,  are 
compared  with  x^p,  1.e.»  >uXn  -  xca]c.  The  mean  of  the  absolute  value  of  xeXn  - 
xca]^  is  Included  in  Table  8-2  for  the  fourteen  sequences.  Also  listed  in  Table  B-2 
are  $he  amounts  of  time  subtracted  from  t  in  Table  A-3  to  obtain  the  time  adjust¬ 
ment  used  in  the  plot  of  Figure  B-2. 

It  is  useful  to  have  an  equation  that  fits  the  shock-wave  versus  x  data  over 
the  range  of  30-150  iren  used  in  Underwater  Sensitivity  Tests.  A  form  that  does 
this  up  to  at  least  120  mm  is 


where 

Rq  =  41.1  mm 
A  =  111.64322 
8  =  1.19226 

The  values  of  A  and  B  are  obtained  from  a  least-squares  fitting  routine,  over  the 
range  of  x  =  10  to  120  nsK,  which  gives  a  QME  of  0.0082  for  the  overall  fit. 

This  QME  indicates  that  the  fit  is  very  good.  A  comparison  of  measured  velocities 
and  predicted  velocities  from  Equation  (B-4)  is  made  in  Table  B~l.  Note  that 
Equation  (B-4),  the  shock  wave  jump  equations,  and  the  equation  of  state  of  water 
are  used,  as  previously  described,  to  give  a  relation  of  peak  shock  pressure  in 
water  as  a  distance  from  the  pentolite  sphere. 

Tt  would  be  useful  to  hydrodynami cs  and  underwater  damage  modelers  to  have  a 
form  for  pressure  versus  scalable  distance,  R/R0.  This  form  would  allow  easy  pre¬ 
dictions  of  shock  propagation  in  water  for  various  size  spheres  of  pentolite.  To 
this  end  a  few  forms  of  shock  velocity  versus  R/Rn  are  evaluated.  Some  additional 
shock-velocity  data  from  pentolite  spheres  of  1.2?  cm  radius  are  included  to  allow 
the  least-squares  fitting  of  these  forms  over  a  range  of  R/R0  of  1.24  to  12.27. 

The  following  equations  are  evaluated. 
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TABLE  B-2.  THE  AVERAGE  DEVIATION  OF  K,  MEAN  RESIDUAL  OF  x, 


THE  TIME  ADJUSTMENT 

FOR  x  = 

0  FOR  FOURTEEN  SEQUENCES 

SEQUENCE 

NO. 

AV  DEV  OF 

If  a 

NO.  OF 
ENTRIES 

1 

i*exp~xcalc' 

1 

0.1588 

11 

0.23 

2 

0.1004 

11 

0.14 

3 

0.1833 

11 

0.27 

4 

0.1229 

13 

0.19 

5 

0.1553 

12 

0.23 

6 

0.1769 

12 

0.26 

7 

0.1129 

10 

0.16 

8 

0.1325 

10 

0.22 

9 

0.1184 

10 

0.19 

10 

0.1687 

9 

0.24 

11 

0.2376 

8 

0.36 

12 

0.2372 

7 

0.37 

13 

0.2700 

5 

0.47 

14 

0.2146 

7 

0.35 

AND 


AT  x  »  0 
(us)b 

-3.7175 
-5.7704 
-3.0573 
-4.7747 
-3.7861 
-5.8653 
-3.4071 
-2! 5869 
-5.0821 
-4.4923 
-3.1908 
-4.3255 
-4.5228 
-3.2966 
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Vco  1+-*- 


us  ■  co  1  + 


A  +  R 


Us  =  c0  I  +  K'  (A  +  l’B 


(A  +  R,)‘ 


where  R1  =  "/R  . 


Table  B-3  gives  the  results  of  least-squares  fitting  these  forms.  The  QME 
for  Equation  (B-5)  is  seen  to  be  large  which  indicates  that  this  form  is  not 
adequate  over  the  entire  range  of  R/Rq,  i.e,9  0  to  «.  However,  it  does  fit  the 
data  very  well  over  the  range  used  in  the  Underwater  Sensitivity  Test.  Equation 
(B-5)  is  implicitly  used  in  the  expression  for  P(R)  in  reference  B-l.  Although 


TABLE  B-3. 


THE  PARAMETERS,  A,  B,  AND  K' ,  FOR  EQUATIONS  (B-5),  (B-6),  AND  (B-7) 


EQUATION 


(B-5) 


2.76523 


1.60598 


1.000000 


0.02934 


(B-6) 


-0.354235 


1.34374 


4.474863 


0.01294 


(B-7) 


-0.360619 


1.28486 


4.435364 


0.01192 


it  does  not  fit  the  data  quite  as  well  as  Equation  (B-4),  it  was  selected 
because,  being  scalable,  it  has  a  much  wider  application.  (Note:  If  one 
derives  A  and  B  from  the  expression  for  P(R)  found  in  reference  B-l,  the  values 
will  be  somewhat  different  from  those  shown  in  Table  B-3.  These  differences 
occur  because  the  fit  is  over  a  shorter  R/R  range  in  the  reference.  In 
addition,  more  recent  data  and  a  more  refined  analysis  technique  are  used  in 
this  report  which  adds  slightly  to  the  differences.) 


Frankel,  M.  J.,  Liddiard,  T.  P.,  and  Forbes,  J.  W.»  "Low-Level  Reaction 
Thresholds  in  High  Explosives  and  Propellants,"  Combustion  and  Flame,  Vol .  45, 
No.  1,  Oan  1982,  p.  35. 
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As  to  the  more  recently  derived  expressions,  the  QME  of  the  fits  for  Equa¬ 
tion  (B-6)  and  (B-7)  are  very  good  and  about  the  same  magnitude  over  the  entire 
range  from  R/R0  =  1  to  ».  This  indicates  that  either  form  could  be  used  to 
reproduce  accurately  the  velocity  over  the  entire  range.  All  of  the  above  equations 
involving  Us  are  forced  to  give  the  sound  velocity  in  water  (1.483  mm/us  at  20°C) 
for  very  large  values  of  Rj.  Equations  (B-6)  and  (B-7)  give  shock  velocities  of 
5.958  and  5.917  mm/us,  respectively,  for  Rj  =  1.  These  extrapolated  values  are 
very  near  the  measured  value  of  6.053  mm/ps. 

It  also  would  be  useful  if  the  shock  velocity  at  Ri  =  1  is  set  equal  to  the 
measured  value  of  6.053  mm/ps.  This  is  easily  done  with  Equations  (B-6)  and  (B-7) 
by  fixing  K' at  3.0816  mm/ps  and  least-squares  fitting  for  A  and  B.  The  results 
of  this  fitting  are  given  in  Table  B-4.  Note  that  the  QME 1 s  of  these  two  param¬ 
eter  fits  are  slightly  higher  than  for  the  three  parameter  fits  given  in  Table  B-3. 


TABLE  B-4. 

THE  PARAMETERS,  A  AND  B,  FOR  EQUATIONS  (B-6)  AND  (B-7)  WITH  K*  =  3.0816 

EQUATION _ A _ B _ QME 

(B-6)  -0.406213  1.30135  0.014Q2 

(B-7)  -0.434278  .1.2231  0.01435 


However,  these  QME's  and  consequential  fits  to  the  data  are  still  quite  good.  To 
obtain  a  better  appreciation  of  these  fits.  Table  B-5  gives  the  measured  data  and 
the  various  calculated  values.  The  two  pararneter  fit  for  Equation  (B-6)  is  chosen 
as  the  best  overall  representation  of  the  data.  This  equation  combined  with 
Equation  (2)  in  the  text  gives  the  peak  shock  pressure  as  a  function  of  R/R0. 
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TABLE  B-5.  A  COMPARISON  OF. FITS  OF  SEVERAL 
OF  THE  EQUATIONS  CONSIDERED 


SHOCK  VELOCITY,  (mm/us) 


R/R0 

MEAS. 

EQ. (B-5) 

EQ. (S-6)a 

EQ. (B-7)a 

EQ.(B-6)b 

EQ.(B-7)b 

1.243 

4.405 

4.375 

4.416 

4.415 

4.429 

4.435 

1.487 

3.649 

3.651 

3.623 

3.626 

3.620 

3.621 

1.730 

3.157 

3.183 

3.149 

3.150 

3.143 

3.141 

1.973 

2.829 

2.860 

2.835 

2.834 

2.830 

2.827 

2.217 

2.596 

2.625 

2.612 

2.611 

2.608 

2.606 

2.460 

2.434 

2.449 

2.447 

2.446 

2.445 

2.443 

2.703 

2.311 

2.313 

2.321 

2.319 

2.320 

2.319 

2.946 

2.218 

2.206 

2.221 

2.220 

2.222 

2.221 

3.433 

2.084 

2.049 

2.074 

2.073 

2.076 

2.077 

3.920 

1.975 

1.940 

1.972 

1.972 

1.975 

1.977 

8.015 

1.686 

1.628 

1.663 

1.666 

1.669 

1.674 

12.269 

1.583 

1.556 

1.584 

1.587. 

1.589 

1.594 

a  LSQ  fitted  parameters  from  Table  B-3 


b  LSQ  fitted  parameters  from  Table  B-4 


APPENDIX  C 


WEAK-SHOCK  SOLUTION  FOR  EXPLOSIVELY  GENERATED  SHOCK  WAVES  IN  WATER 


C_i 

In  the  work  of  Rogers  it  is  shown  that  the  peak  pressure  of  a  weak  shock 
wave  is  given  by  the  expression 


(Cl  +  2(R,/jl  )  InCR/R,)]*5  -  1} 
P  =  P1  Wa0)  MR/Rjj 


(C-l) 


where 

Pj  =  the  peak  pressure  of  the  shock  at  distance,  Rp  from  center  of  charge 

3 

l  =  the  characteristic  length  =  p.  c^  t./P, 6 
0  0  0  1  i- 

pQ  =  the  initial  density 

cQ  =  the  sound  velocity 

3  =  the  parameter  of  non-linearity  of  a  fluid. 

The  expression  for  the  pulse  time  constant  is 

x  =  Tj  [1  +  2(R1/Aq)  ln(R/R1)DJs  (C-2) 

It  is  not  obvious  that  these  two  equations  scale  for  explosive  charges  of 
different  initial  sizes.  However,  the  gage  data  of  Table  4  in  the  text  show  that 
Pj  and  t/Rq  are  approximately  constant  at  the  scaled  distance,  R l / R0 .  For  a 
reference  shock  pulse  (Pj,  Ri,  q),  Ri/Rq  also  is  a  constant.  Then,  £0  -  tjKq  = 
R0Ki  =  RjKg,  where  K0  *  pc03/Pi,  Kj  and  K?  are  unknown  constants  for  a  given 
reference  state  (Pj,  Rj,-q).  Substituting  these  simplifications  into  Equation 
(C-l)  results  in 


P  =  {[1  +  2K3  ln(R/R1) -  1} 
T  (R/R1)'K3  ln(R7Rp 


(C-3) 


^Rogers,  P.  H.,  "Weak-Shock  Solution  for  Underwater  Explosive  Shock  Waves," 
Journal  of  the  Acoustical  Society  of  America,  Vol.  62,  1977,  pp.  1412-1419. 
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x  =  x1ll  +  2K3  ln(R/R1)3  15  (C-4) 

where  K3  =  1/K„ ,  an  unknown  factor.  Equations  (C-3)  and  (C-4)  are  preferable 
when  fitting  to  weak-shock  data.  Note  that  these  equations  were  derived  using 
weak-shock  approximations  and  the  assumption  that  the  pulse  shape  is  exponentially 
decaying  in  time.  These  assumptions  limit  the  general  usefulness  of  the  equations, 
but,  the  work  of  Rogers  is  a  large  step  forward  in  treating  weak  shocks  in  water. 

Tables  C-l  and  C-2  give  the  results  of  fitting  the  shock-velocity  data  of 
Table  2  in  the  text  to  Equation  (C-3)  for  two  different  reference  shock  pulses.  It 
is  clear  that  Equation  (C-3)  gives  a  fair  fit  to  the  peak  pressure  data  for 
Pj  =  19.1  and  14.1  kbar.  Equation  (C-3)  fits  the  data  better  for  the  lower  refer¬ 
ence  shock  with  peak  pressure  of  14.1  kbar.  It  is  expected  that  Equation  (C-3) 
will  closely  reproduce  the  peak  pressure  data  for  pressures  less  than  a  few  kilo- 
bars. 

Equation  (C-3)  also  was  fitted  to  the  lower  pressure  gage  data  of  Table  4  in 
the  text.  The  results  of  a  least-squares  fit  of  Equation  (C-3)  to  the  data  is 
given  in  Table  C-3.  As  expected,  this  fit  is  better  tha,  .hat  indicated  by  the 
results  given  in  Table  C-l  or  Table  C-2,  since  it  was  maoc  over  a  lower  pressure 
range. 

The  decay  time  (pulse  half-width),  x,  given  in  Table  4,  is  the  time  taken  for 
the  pressure  to  fall  one-half  of  its  peak  value.  Equation  (C-4)  can  be  used  to 
calculate  this  x  if  the  shock  decay  is  exponential  and  xQ  is  chosen  as  the  pulse 
width  at  half  the  peak  pressure.  Note  that  the  constant,  K3  =  1.9,  was  obtained 
by  fitting  the  pressure  data.  The  calculated  values  of  x/R0  are  given  in 
Table  C-3,  along  with  the  measured  x/R0  data  from  Table  4.  These  values  are 
within  10  to  20%  of  the  measured  values.  An  improvement  in  reproducing  the  data, 
using  Equation  (C-4),  would  be  expected  if  more  realistic  shock  profiles  were  used 
in  the  derivation  of  Equation  (C-4). 


_o _ (kbar)* _  (kbar) _ (kbar) 


2.217 

14.3 

13.89 

0.41. 

2.460 

11.4 

11.13 

0.27 

2.703 

9.4 

9.33 

0.07 

2.946 

8.0 

8.04 

-0.04 

3.433 

6.1 

6,29 

-0.19 

3.920 

4.7 

5.15 

-0.45 

8.02 

1.63 

1.97 

-0.34 

12.27 

0.75 

1.17 

-0.42 

*0btained  from  Table  2  in 

text. 

TABLE  C-2.  WEAK  SHOCK-WAVE  PEAK  PRESSURE  FOR 
P1  =  11.4,  R1/RQ  =  2.46,  AND  K3  =  2.9 

R 

Ro 

P 

(kbar)* 

Pcalc 

(kbar) 

P~Pca1c 

(kbar) 

2.703 

9.4 

9.25 

0.15 

2.946 

8.0 

7.83 

0.17 

3.433 

6.1 

6.02 

0.08 

3.920 

4.7 

4.89 

-0.19 

8.02 

1.63 

1.84 

-0.21 

12.27 

0.75 

Obtained  from  Table  2 

1.09 

-0.34 

C-3 
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TABLE  C-3.  MEASURED  SHOCK-WAVE  PROFILES  IN  WATER 
FROM  THREE  SIZES  OF  PENTOLITE  SPHERES 


R  MEAS.  PK,  CALC.  PK.  MEAS.  CALC. 

1L  PRESSURE  PRESSURE  t/R0  t/R0 

^  (kbar)  (kbar)*  (us/cm)  (us/cm) 


4.97 

4.11 

3.048 

3.052 

6.00 

7.25 

6.50 

6.84 

2.038 

2.003 

8.36 

9.44 

8.01 

4.11 

1.543 

1.490 

11.73 

10.84 

6.84 

1.449 

1.490 

3.63 

10.84 

1.27 

1.447 

1.487 

9.83 

10.85 

12.27 

1.27 

0.803 

0.851 

12.01 

13.25 

12.28 

4.11 

0.817 

0.850 

12.18 

13.25 

^Reference 
and  K3  =  1 

state  of  Pi 
.9 

=  3.823  kbar. 

R]/R0  =  4.39, 

lt/R0  =  ! 
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